NLO QCD corrections to top quark pair production in association 
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Abstract 



Oh. 

' We compute the QCD corrections to the production of a tt pair in association with one hard jet at 

the Tevatron and the LHC, using the method of generahzed Z3-dimensional unitarity. Top quark 
I decays are included at leading order in perturbative QCD. We present kinematic distributions of 

top quark decay products in lepton plus jets and dilepton final states at the Tevatron and the LHC, 
using realistic selection cuts. We confirm a strong reduction of the top quark forward-backward 
^ ^. asymmetry in pp — > tij at next-to-leading order, first observed by Dittmaier, Uwer and Weinzierl. 

We argue that there is a natural way to understand this reduction and that it does not imply a 
breakdown of the perturbative expansion for the asymmetry. 



1. Introduction 



(N 
> 

j. ■ Experiments at the Large Hadron Collider (LHC) are in the process of testing the very suc- 

I cessful paradigm of the past thirty years - the Standard Model of particle physics - in an entirely 

O^l ' new energy regime. The goal of these experiments is to unravel the mechanism of electroweak 

I symmetry breaking. Electroweak symmetry breaking may occur due to the Higgs mechanism, as 

' in the Standard Model, or it may be a consequence of more complex dynamics. In the Standard 

I Model the unusually large mass of the top quark is introduced by hand and top quarks do not 

play any role in electroweak symmetry breaking. This may change in extensions of the Standard 
Model. Moreover, the large mass of the top quark makes its interactions with an agent of elec- 
troweak symmetry breaking parametrically enhanced, relative to other quarks. Therefore, top 
quarks provide a window of opportunity to investigate physics of electroweak symmetry breaking; 
; ^ ' this motivates a rigorous experimental program to study top quarks and their properties. 

\ More than a decade after the discovery of the top quark in 1995 [1, iS], the Tevatron remains at 

the forefront of top quark studies. During that time, CDF and DO reached remarkable precision 
in measuring the top quark production cross-section in pp collisions and the top quark mass Q . 
Also, the top quark charge is constrained by experimental measurements i^|5|] and the hypothesis 
that the top quark spin is one half follows indirectly from the pp — ^ tt cross-section measurement 
and is consistent with the observation of the top quark decay t — > Wb. The structure and the 
strength of tbW interaction vertex was thoroughly studied [6, 7] and limits on exotic contributions 
to top quark decays were placed. A recent measurement of the forward-backward asymmetry of 
the top quark 8, 9] may point towards non-standard contributions to top quark pair production; 
currently, such a possibility is being actively discussed. 

In spite of the enormous integrated luminosity accumulated by the Tevatron, many of their 
measurements are limited by statistics since the ratio of the hadron center of mass collision energy 
to twice the top quark mass is not large. This will change with the start-up of the LHC since the 
LHC will be a tt factory. It is expected that many of the studies performed at the Tevatron will 
be continued and extended at the LHC with much higher statistics. The higher energy of the LHC 
leads to a much larger cross-section for the top quark pair production, compared to the Tevatron, 
and, also, to a greater variety of the kinematics of the produced top quarks. In contrast to the 
Tevatron, top quarks are produced with larger energies and transverse momenta at the LHC. This 
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kinematic feature implies larger probability for top quarks to radiate gluons and the increase of 
the relative importance of the + jets final state in the inclusive tt sample. 

It is important to understand how many jets are produced in association with a tt pair. Indeed, 
since top quarks may decay into leptons, missing energy and jets, tt production is an important 
background to standard signatures of physics beyond the Standard Model, where the number 
of jets in the process is often taken as a discriminator against the background. Note also that 
selection cuts that enhance New Physics contribution to a particular observable typically select 
harder top quarks, thereby enhancing the probability of QCD radiation. An interesting example is 
the study of background to the Higgs boson production in vector boson fusion pp H+2 jets — > 
W^W~ + 2 jets — >■ lliyD -\- 2 jets. The final state is identical to the ti production if both top quarks 
decay leptonically and the two jets are associated with 6-jets from top quark decays. Requirement 
that the two jets are separated by a large rapidity interval suppresses the pp — >■ tt background 
since 6-jets from top quark decays tend to be central. On the other hand, in pp — > + jet process 
there exists kinematic configuration with one central, one forward and one backward jet, separated 
by large rapidity gap. Since one of the forward/backward jets can be a gluon jet from the initial 
state radiation, only one 6-jet from top decays has to be non-central. Because the probability to 
have one &-jet from top decays going in the forward/backward direction is much higher than the 
probability that 6-jets from both top decays are in the forward/backward direction, tt+jet becomes 
the dominant background. There are other cases in which ti + jets background is important and 
tt+jet is an essential part of it. For example, the importance of the tt+jets background increases 
with the number of jets [IQj . for supersymmetry searches in events with jets + missing energy. 

Accurate theoretical description of pp(pp) — )■ tt + jets process requires next-to-leading order 
(NLO) QCD calculations. Such calculations are difficult because of the large number of diagrams 
and the presence of massive particles. While NLO QCD corrections to ti production are known 
since early 1990s 11, 12, l^[l3|, a calculation of NLO QCD corrections to ti -I- jet was reported 
recently 3; 16 1. Even more recently first results for NLO QCD corrections to pp — ti + 2 jets 
were presented in Ref. NLO QCD corrections are also known for other associated production 



processes that involve top quarks, including tij [18|, tiZ [l9|, tiH [20| and tibb [2l|, |22|, |25 
However, quite often, when NLO QCD corrections are computed, top quarks are treated as stable 
particles and their decays, even at leading order, are not included. This is not entirely realistic 
since all cuts designed to either select a ti sample or to reject it, to search for physics beyond 
the Standard Model, apply to top quark decay products. It is well-known that kinematics of the 
decay products is affected by spin correlations. Therefore, it may be important to account for 
decays of top quarks in the calculations of NLO QCD corrections, to have full confidence in the 
results. We note in this regard that all acceptances employed in experimental studies of top quarks 
convert a measurement with cuts on the top quark decay products to a result that refers to top 
quarks as stable particles, are currently calculated with parton shower event generators and are 
never corrected for higher order QCD effects. This is equivalent to declaring, without any proof 
or argument, that QCD radiative corrections to top quark pair production are always reduced to 
kinematic-independent i^T-factors and, therefore, cancel out in the calculation of all acceptances. 
Whether or not such an assertion is correct, can only be established through a next-to-leading 
order calculation that includes top quark decay products. We note in this regard that computation 
of NLO QCD corrections to observables in ti pair production that are sensitive to top quark decay 
products was pioneered in Refs. 24, H, 0, [i^, S, [i^]. The complete calculation of NLO QCD 



corrections to ti pair production that retains spin correlations, allows application of arbitrary cuts 
to the decay products of top quarks and includes QCD corrections to the decay, was presented 
very recently in Refs. 3^ 31 1. 



In this paper, we extend the results of Refs. [l^ 13 by computing NLO QCD corrections 
to pp(pp) ti + ]et, including leptonic and hadronic decays of top quarks at leading order. 
To calculate one-loop virtual amplitudes, we employ the method of generalized i?-dimensional 
unitarity suggested in Ref. [sl] and applied to massive particles in Ref. [s^. We point out that a 
basic object that needs to be calculated for each phase-space point within the unitarity framework 
is a one-loop helicity amplitude. Since helicity amplitudes need to be calculated anyhow, it takes 
(roughly) the same effort to calculate helicity density matrix p\x> ^ ^\^x' i needed for an exact 
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account of top quark decays, and the spin-averaged cross-section a ^ ^\P\x- Therefore, when 
unitarity methods are employed in the calculation of one-loop corrections, there is no good reason 
to avoid including the leading order on-shell top quark decay. Note, however, that in order to 
produce a result for pp{pp) — > tt -f jet where NLO QCD corrections to both production and decay 
are included, one needs to know NLO QCD corrections to the decay t — > bW +iet. Such calculation 
is straightforward but tedious; we hope to return to it in a separate publication. 

The remainder of the paper is organized as follows. In Section [2] we comment on some aspects 
of generalized unitarity that are important for the calculation that we report in this paper. The 
goal is to provide details which were not presented clearly in the original publications [32, [s^ . In 
Section [3] technical aspects of the calculation are reviewed. In Section |4] phenomenological results 
are described. We conclude in Section [SJ 

2. Aspects of generalized D-dimensional unitarity 

To calculate one-loop virtual corrections, we employ the framework of generalized D-dimensional 



unitarity, following Refs. [32|, [33|. This method grew out of the observation of Ref. [3^ that 



the procedure for tensor reduction of one- loop integrals, suggested in Ref. [35|, meshes well with 



unitarity-based methods and allows computation of one-loop on-shell scattering amplitudes, rather 



than Feynman diagrams. The central idea of Refs. [3^ |33[ is that this technique can be extended 
in such a way that both, cut-constructible and rational parts of the amplitude, can be obtained 
in an efficient way. 

While Refs. (sil . [33j explain clearly how numerical computation can be set up, it is beneficial 
to discuss some aspects of the implementation in more detail. Recall that within the context of 
generalized I?-dimensional unitarity, we deal with the so-called primitive amplitudes [36] . which 
are gauge-invariant subsets of color-ordered amplitudes. The importance of primitive amplitudes 
for generalized unitarity follows from the fact that external particles in a primitive amplitude are 
ordered and no permutations are allowed. This feature enables enumeration of all the propagators, 
that contribute to such an amplitude, in an unique fashiorQ. 



As explained in Refs. |32l. I33|. a full one- loop amplitude can be computed if the integrand of 
the one-loop amplitude is known in the extension of the four-dimensional theory to 13 > 4 integer 
dimensions, and the loop momentum is restricted to five dimensions. In this case, an integrand 
of any iV-point primitive amplitude can be decomposed into a sum of terms with at most five 
Feynman propagators 

fj (j c b n 

A n\ \ ^ iii2i3i4,h i \ ^ "'111213*4 , \ ^ ^111213 , \ ^ ili2 , \ ^ il /1 \ 

- d.,d.MA. ^ dMdZ + ^ V, + d-d- + rf-- W 

[hKs] [nyii] [»i|»3] [n|22] [n] 

In Eq.([T]) we use the notation [ii|i„] = 1 < ii < 12 <...<«„< iV and di{l) = {I + PiY — mf, 
i = 1...N. The left hand side of Eq.(IT]) is completely specified by Feynman rules. The challenge is 
to find an efficient and numerically stable way to find the coefficients e^^ ^5 j on the right hand 
side of Eq. ([T]) . It was suggested in Ref. [35j that a very efficient way to do that is to solve for the 
coefficients e,^ j^,, a^^ sequentially, by calculating left and right hand sides of Eq.([T]) for special 
values of the loop momentum I, where certain subsets of inverse Feynman propagators di, c?2, d]\[ 



vanish. All propagator sets with up to five members should be considered. It was shown in Ref. [3J] , 
that the coefficients e^^ ir j have restricted functional dependence on the loop momentum I; 

they are typically given by linear combinations of constant terms and traceless tensors of ranks 
up to five, four, three, two and one, respectively. Those tensors are defined on linear spaces 
constructed from basis vectors, orthogonal to momenta present in Feynman propagators of a 
corresponding term in Eq.([T]). The dimensionality of these transverse spaces is — 4, D — 3, 
D — 2, D — 1, D for terms with five, four, three, two and one denominators, respectively. It can be 



^Recently, a procedure that avoids the necessity to deal with primitive amphtudes in the context of unitarity- 
based one-loop computations was described in Ref. [37| . 
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shown that, after integration over the loop momentum, the contribution of the traceless tensors 
vanishes. 

To illustrate this, we consider a term with two denominators di-^^i^ = di^2 in Eq.(IT|). We choose 
the momentum I in such a way that di_2 read 

di=f-ml, d2 = [l + pf - ml. (2) 

We write I as a linear combination of the vector p, that enters d2 and the vector l±_ that belongs 
to the transverse space described earlier 

It" = xp" + li^; pl^^O. (3) 

The numerator function 612(0 given by the sum of two terms - an /-independent constant 6^2^ 
and a function 612 that depends on the transverse component of the momentum / 

bi2{l)^bf^^ +h2{l±). (4) 
The integral of the function &i2(/j_) over directions of the transverse space vanishes 



d''-'l^5{li-^It,)bl2(l±)^0■ (5) 
Because di,2 depend on ll_, Eq.® implies that bi2{l±) does not contribute to the final result where 



one- loop amplitude is expressed through scalar (master) integrals. As explained in Refs. [35, 

|32, 33], renormalizability of the theory and the fact that momenta of all external particles are kept 
in four-dimensions restricts the functions eij...^^, . For example, as we already mentioned, 
ranks of tensors that contribute to those functions are limited and, also, those functions must be 
parity even with respect to the fifth component of the loop momentum ^5 . 

The general set up that we just presented is employed in many existing applications of general- 
ized unitarity is discussed in detail in the original papers [32I, fssj. However, 
there are some aspects of the implementation of this method that are not adequately described in 
those references and that warrant some discussion. This is what we do in the remainder of this 
Section. 

Our first comment concerns the coefficient of the term with five denominators ei-^i^.A^- Accord- 
ing to the preceding discussion, it can be any even tensor up to rank five, composed of I5 . Hence, 
the most general form of this function is 

e = eo -(- ei?5 -I- 62/5. (6) 

However, it is easy to see that 

Il^fo + O{di,d2,...d5), It^ fS + 0{di,d2,...d5), (7) 

where /o is an /-independent term that only depends on the kinematic variables of a particular 
five-point function. It follows from Eq.© that either scalar five-point integral or the five-point 
integrals with /g or /| in the numerator can be chosen to be the five-point master integral since 
the difference between these integrals is given by linear combinations of the four-point functions. 
In R ef. |32l| the scalar five-point integral was chosen as the master integral but, as was realized 
later [38|, that choice is somewhat unfortunate. Indeed, it is well-known that if one calculates the 
cut-constructible part of a one-loop amplitude, one does not need to consider the five-point master 
integrals. On the other hand, if one chooses e to be an /-independent constant, the scalar five-point 
function appears as a master integral. Because in the I? — > 4 limit the five-point function is not 
independent of the four-point functions, and because this is the limit of interest, large numerical 
cancellations between the four-point and the five-point master integrals are to be expected and 
indeed occur. These cancellations may be so strong that the accuracy of the result deteriorates. 
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To ensure that all numerical cancellations happen locally, it is useful to choose the basis of master 
integrals such that ei^..^^ ~ ^| ~ 'i)-4- Since 

f d^l 1% 4 

the new master integral does not contribute to the final result, but it is needed at the intermediate 
stages for proper identification of lower-point functions through the residues of the one-loop scat- 
tering amplitude. The new choice of the master integral for the five-point function helps to avoid 
large numerical cancellations. A similar discussion of the role of the five-point master integral in 
the context of Z?-dimensional unitarity was recently given in Ref. [40| . 

Next, we describe a procedure to find functions ei^..^^, ...,aij in a numerically stable way. We 
put special emphasis on exceptional cases where the general reduction algorithm needs to be 
extended. To simplify the discussion, we consider four dimensional unitarity and focus on terms 
with two denominators in Eq.([T]). We stress that all the subtleties associated with the reduction 
can be illustrated by those considerations. To this end, we pick any term with two denominators 
in Eq.([T|); we refer to the inverse Feynman propagators in that term as di and d2 and to its 
numerator as For ease of the presentation, we only discuss the cut-constructible part in what 
follows. We assume that di and d2 and the momentum / are given by Eqs. (|2|3|) and write l±_ as 



/_L = ^ XiVi, pVi = 0, ViVj = 6ij. (9) 

In Eq.®, we introduced basis vectors Vi, i = 2, 3, 4 for the linear space which is transverse to the 
momentum p. As explained in (s^. [ssj. the function b{l) can be written as 

b{l) = &i + b2{v2l) + bsiv^l) + b^ivil) + 65 {{V2lf - (Vslf) 

+b6 {{V2lf + {V3I? ~ 2{vAlf) + bj{v2l){v3l) + bs{v2l){v4l) + bg{v3l){vAl) . (10) 

In Ea. (jlOp bi g are unknown ^-independent parameters; the goal is to set up an algorithm to find 
them. To calculate foi,.,9, we would like to calculate b{l) for some values of the loop momentum 
and then use the results of the computation to solve the system of linear equations where 61. ..9 
are treated as the unknown parameters. 

It was suggested in Ref. that it is convenient to choose the momentum / for which both 
inverse Feynman propagators vanish di{l) = d2{l) — 0. In Ref. [s^l it was pointed out that, for 
such choices of the loop momenta, the function b{l) and all other numerator functions in Eq.(IT]) 
can be calculated using tree-level on-shell scattering amplitudes. This feature is apparent from 
Eq. ([ij since functions e, d, c, 6, a are the residues of one-loop scattering amplitudes in a situation 
where certain virtual particles go on the mass-shell. Such residues are computed by cutting some 
internal lines in an one-loop amplitude; this terminates the flow of the loop momentum and turns 
a one-loop amplitude into a product of tree-level on-shell amplitudes. These on-shell amplitudes 
are conventional apart from the fact that they must be calculated for the external complex on-shell 
momentum. Using Eq.(j9l), we find that components of the momentum I for which di^2 = are 
subject to the following constraints 

{ml - ml - p'^) ,2222 

V ' 'j- = "ii-2;p- (11) 

It follows from Eq. ([TT|) that there are infinitely many momenta I that satisfy the di.2 = condition 
since only l']_ is constrained. Therefore, all we need to do is to choose nine vectors with different 
directions of l±, calculate b{l) for those vectors and solve the resulting linear system of equations 
to find &i,.,9. 

In principle, the nine vectors can be chosen arbitrarily. However, a more systematic procedure 
emerges if we parameterize the transverse component of the loop momentum as 

l± — \l± \ (sin 9 cos (j) V2 + sin 9 sin (p V3 + cos 9 V4) ■ (12) 
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When this parameterization is employed in Eg. dTUl) . b{l) becomes a degree- four polynomial of 
e'**, e*"^. Since we are interested in finding coefficients of that polynomial, the technique of discrete 
Fourier transform can be employed. In the context of unitarity methods, this technique was 



described in Refs. [42|, |43| and we do not repeat it here. 

Note, however, that the application of the discrete Fourier transform requires division by 
According to Eq. ([TT|) . vanishes if = x^p^ which corresponds to = (to2 — toi)^ or 
p^ — (to2 + mi)'^. These kinematic points are not dangerous if only massless virtual particles 
are considered since the offending integrals are scale-less two-point functions that are discarded 
in dimensional regularization. The situation changes once virtual massive particles are considered 
and this is the case that we are interested in. Note also that close to those exceptional values of 
p^, can be small, so that division by it may lead to numerical instabilities. 

To handle the case of small \l±\ in a numerically stable way, the method of discrete Fourier 
transform is not directly applicable and the system of equations must be solved differently. There 
are many ways to solve a system of linear equations avoiding division by the procedure 
that we have implemented in the numerical program is described below. We begin by choosing 

= x±V2 i ■ = l^-Lp- Recall that l'^ is fixed from the on-shell condition Eg. dTTI) and 

therefore 2:3 is expressed through x±, X3 = \/\l±\'^ — a;5_. We calculate b± = b{l'^) and remove x§ 
in favor of and x^ where possible. We obtain 

b± = bi + b2Xj_ ± X3b3 + 65 {2x1 - ll) + bell ± b7X±X3. (13) 
Taking the sum and the difference of b± , we arrive at 

^^±±M = b, + b2X^ + b, {2x1 -ll)+b,ll, ^^^^^^h3 + hx^. (14) 

The right hand sides of these equations are polynomials in Therefore, we can apply discrete 
Fourier transform with respect to x±^ to find coefficients 62, ^5, b"^ = b\ + b%l\ as well as 63, 67 in 
Eq.dHl). 

To determine the remaining coefficients, we take the vector Z_l to be in the V2 — V4 plane. We 
choos^ X2 = y^l'^~ I and consider three different vectors Z^-* = X2V2 + W4, Z^"* — —X2V2 + Vi, 

= — a;2W2 ^ We use the notation ba — b{xp + Z^'') but we assume that all terms with 
previously computed coefficients 
that the following relations hold 



previously computed coefficients are subtracted when b(xp + /^^) is computed. It is easy to see 



{ba-b,) ,_(&.-M {b,-bf-b4-bsVK^ 



(15) 



Also, we find bi = bf^ — l^b^ and, finally, the term with bg can be obtained by calculating the 
function b for the momentum l± ~ V3 + ^y\l±^^-l Vi. 

We have just described a method to calculate coefficients 61. .9 in a numerically stable way for 
arbitrary values of |Zj^|. In the numerical program, we switch from the discrete Fourier transform 
to the solution just described, depending on the value of \l±\. However, the described methods can 
only work if the decomposition of the loop momentum, as in Eq.Q, exists. A glance at Eq. (|lip 
makes it clear that the decomposition fails for the light-like momentum, p^ = 0, and we have to 
handle this case differently. We describe a possible solution below. 

Because we are interested in one-loop calculations for infra-red safe observables, it is reasonable 
to assume that the vector p can be exactly light-like but it is impossible for that vector to be nearly 
light-like, since such kinematic configurations are, typically, rejected by cut^. Hence, we have to 



^One should not be confused by the mismatch in dimensions. In an actual computation, we always rescale all 
dimension-full variables by the center of mass collision energy, to turn them into order one numbers. 

^External particles with small masses are obvious exceptions but rarely do we need to know observables for, say, 
massive b-quarks in a situation when all kinematic invariants are large. 
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modify the above analysis to allow for a light-like external vector. To this end, we choose a frame 
where the four-vector in Eq.(I2]) reads p = {E,Q,Q,E). We introduce a complimentary light-like 
vector p = {E, 0, 0, —E). The loop momentum is parameterized as / = xip + X2P -f We denote 
the basis vectors of the transverse space as W3,4; they satisfy viVj — 5ij, PW3.4 = 0, pf3.4 = 0. The 
on-shell condition for the loop momentum fixes X2 



2 2 

TOn — TOT „ , 

X2 = — s = 2pp, (16) 



and a linear combination of Xi and 

+ mlxi ~ ml{l + xi) = 0. (17) 

The parameterization of the function b reads 

b{l) - 61 + 62(^0 + b3{v3l) + b^iv^l) + b^{pl){pl) + be{pl){vsl) 

+ bj{pl){Vil) + bsiM^ - (v^lf) + bg{v3l){v4l). (18) 

We describe a procedure to find the coefficients 5i, ..,&g in a numerically stable way. To this 
end, we choose xi = 0.5. This fixes and X2 is fixed by the on-shell condition Eg. dTBl) . The 
freedom remains to choose the direction of the vector l±_ in the V3 — W4 plane. Consider four 
different vectors 

z'l ^ = + X4V4, 1^^ — -~V3 + X4U4, I'f^ = V3 — X4V4, if^ = — ■i;3 — a;4t;4, (19) 

where X4 — — 1. We use vectors /^"^ = x\p^X2P + l^'^ ■, a — a,b, c, d, to calculate the function 
= 6(P). Using ba,.:b d, we can immediately find the coefficient 69 

For the determination of the remaining coefficients, it is convenient to introduce two linear com- 
binations 

647 = {b^^'' + - b^'^) ' bse^l {b^"'' - M")) - &9v/zi~I. (21) 

As the next step, we choose xi = —0.5. Note that this changes the value of according 
to Eg. dTTl) . We then repeat the calculation described above. Our choices of momenta in the 
transverse plane l± are the same as in Ea. (|19p but, to avoid confusion, we emphasize that X4 has 
to be calculated with the new We will refer to b computed with those new vectors as b^'^\ 
b'^^\ etc. We calculate 647,36 by substituting b^"^ — > 5" in Eq. (PT|l . It is easy to see that simple 
linear combinations give the desired coefficients 

64 = i (b47 + ^47) , ^7 ^ (b47 - hi 

bs^l (b36 + he) , 65 = - (bse ^he] ■ (22) 



V / s 

Other coefficients, required for the complete parameterization of the function b{l) in Ea. ([T8)) . are 
obtained along similar lines; we do not discuss this further. However, we emphasize that the 
procedure that we just described is important for the computation of one-loop virtual amplitudes 
in a situation where both massless and massive particles are involved. In particular, it is heavily 
used in the computation of NLO QCD corrections to top quark pair production discussed in 
Ref. [i^l and in this paper. 
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As a final remark, we note that there is another important difference between reducing the 
two-point function to scalar integrals for a light-like and a regular vector p. Indeed, the two-point 
function with 7^ can be immediately simplified 



(27r)^did2 J (27r)^did2 

This feature follows from the fact that all but the /-independent terms in the function b{l), Ea. pU)) . 
vanish when integrated over the directions of the transverse space, cf. Eq.([S]). Hence, the only 
integral we need to know is the scalar two-point function. However, in case of a light- like vector 
= 0, three master integrals survive even after averaging over the directions of the vector I in 
the (two-dimensional) transverse space 

d^i b{i) _ f d^i b, + b2{pi) + b5{piy_ ^24) 



{2Tr)D did2 J {2tt)D did2 

Those integrals must be included to the basis of master integrals in case when double cuts with 
a light-like vector are considered. The calculation of those additional master integrals is straight- 
forward; for completeness, we give the results below for the equal mass case mi — ni2 — m. We 
introduce di — {I + fci)^ — m^, d2 = {I + fc2)^ — m^, ^2 — fci = p, p^ — 0, pki — ri, pk2 — ^2, 
cr = {AT:y-^T{l + e)r(l - e)^/T{l - 2e) and find (D = 4 - 2e) 

f d^l Ip _ (ri + rs) / 1 



+ h (25) 



icr J (27r)^ did2 2 \ e \ fi^ 

fi'^ f d^l mm {rl+nr2+4} ( 1 .... , , 



3. Technical details of the calculation 



In this Section, we collect all the information pertinent to the technical details of the com- 
putation. As with any next-to-leading order calculation, we need to deal with one-loop virtual 
corrections, real-emission corrections and the soft-coUinear subtraction terms. Each of these con- 
tributions requires a different treatment and we discuss them in turn. Much of what we say in 
this Section can be found in the literature but we decided to summarize all the required details in 
this paper, for completeness. 

3.1. The decay of the top quark 

As we pointed out in the Introduction, we include decays of the top quark, albeit in leading 
order in perturbative QCD. The way this is done can be illustrated by considering tree- level 
process; NLO QCD effects do not introduce additional subtleties. Our discussion follows Ref. (30| . 
Consider production of a tt pair and a jet in the collision of partons a and 6. Top quarks decay 
into bottom quarks and VF-bosons. For definiteness we consider leptonic decays of the W bosons 
ab — ttj — > {bl~v){bl^i')i but there is not much difference with hadronic decays as long as we 
restrict ourselves to leading order in top decays. The amplitude for this process is written as 

A = ^a{t* bl+v) Aai3{ab t*i* + j) Ep{i* -> bl-p), (27) 

where ^! and S are the off-shell fermion currents and a, /3 are the Dirac indices. The fermion 
currents contain t(t) propagators; for example 

^ bl+iy) = A^{t* ^ bl+i,) \ , (28) 

where A^ is a sub-amplitude that describes transition of an off-shell fermion t* to a final state 
bl+v. 
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We note that radiation of a jet off the top quark decay products is not included in Eq. ipT)) . 
Therefore, as written, Ea. (|27p is not gauge invariant. However, this problem disappears once 
we take the on-sheh limit for both t and i. To take the on-shell limit, we imagine squaring 
the amplitude A and integrating over the phase-space of the final state particles. In the limit 
Tt/mt — > 0, propagators that appear in Eq. (|28l27p become proportional to the delta- function 



1 



2mtr. 



5(p?-m?). (29) 



Finally, we factorize the phase-space into a phase-space for ttj and a phase-space for top decay 
products and observe that the (5- function in Eg . ((291 forces top quarks on their mass-shells. This 
separates production and decay stages and ensures that those stages are separately gauge invariant. 
The top decays are then conveniently implemented by making the substitution in Eq. ()27p 



E^{t^ bl-D) ^ V^ipi) = lt^^±:^^A^{t~^ hl'v). (30) 

V 2mtl t 

We emphasize that U and V are conventional spinors with definite polarizations that are fully 
determined by momenta and helicities of the top quark decay products. In fact, direct multipli- 
cation of spinors and Dirac matrices is an efficient way to reconstruct U and in a numerical 
program for each phase-space point. Finally, we note that similar considerations work for one- loop 
corrections and the subtraction terms. 

3.2. One-loop amplitudes 

Our calculation of the one-loop amplitudes is based on the method of generalized Z3-dimensional 
unitarity. In this Section, we describe the building blocks that enter the calculation. 

The one-loop amplitudes, that we compute in this paper, require renormalization of the bare 
top quark mass rrit barc — Zmmt , the bare top quark field V't.barc = ^f^ipt and the strong coupling 
constant gs. For massive quarks, the mass renormalization and the wave function renormalization 
are performed on-shell. The renormalization constants depend on the regularization scheme; in 
this paper all the results for one-loop amplitudes are given in the four-dimensional helicity (FDH) 



scheme, introduced in Refs. [44|, |45|. In that scheme, the on-shell renormalization constants read 

= ^2 = 1 - Cps'cr + 5^+O(5^e), (31) 

where Cp = [N'^ — l)/{2Nf,), N^. = 3, is the conventional QCD color factor. Finally, we note that 
the renormalization of the coupling constant is the same as in the MS-scheme but a finite shift 
needs to be applied [i^ to relate and conventional 

a, ^ + ——y (32) 

Objects that we compute within the unitarity framework are the so-called primitive amplitudes 
(sij . Regular amplitudes and primitive amplitudes are related through the color decomposition 
procedure. Below we give the decomposition of all amplitudes required for our computation. There 
are two types of tree- and one-loop amplitudes that need to be considered to calculate production 
of a top quark pair and a jet in hadron collisions: tt + ggg and — > tt + qq + g. We first 
consider the two-quark (n — 2)-gluon process. At tree- level, the — > t + t + {n—2) ghions scattering 
amplitude reads 

■AT = E (^"'''' -^"''"^ t ih, 2u a(3)„ a{n)g) , (33) 
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where A^^^^ is the tree (left) primitive amphtude and Sn-2 is the permutation group of (n — 2) 
elements. The SU(3) generators are normalized as Tt{T°'T^) ~ S"-^ and satisfy the commutation 
relations 

[T'',T'']^ ^F^^T". (34) 

This normalization allows us to use color-stripped Feynman rules 47 , 4^ 4^ to compute primitive 
amplitudes. 

A similar expression for the one-loop amplitude is more complicated. Using the color basis of 
Ref. (soj l , we can write the one- loop amplitude as a linear combination of left primitive amplitudes 



^i-ioop _ ^2 (r'^^r^'a ...T^'pT^^)*^ (F""p+i ...F""" ) 

p=2 aeS„_2 

X (-l)"^' (1,, a{p)g, a{3)g,2t, a{n)g, <j{p + 1),) 

1 ""^ / 
+ ^ Grgf (a3,...,a„) n;4/^-^l(l,-,2,,a(3),,..,a(n),) 

+4Y/'*'(lt.2t,a(3)„...,a(n),)Y 



(35) 



where nf = 5 is the number of quark flavors that are considered massless. In Eq. psp . for p — 2, the 
factor {T...T)'^ (T-^^x^i)^ and for p = n the factor {F...F)^^^^ b^^^^. In the second term in 



[1/2,9] 



Eq.([55|). Sn-^j is the subgroup of S'„_2 which leaves GrJ^''^-^ invariant. Primitive amplitudes 

^L-j^'*' contain one closed loop of massless and massive (top quarks) fermions, respectively. The 
color factors read 

Grgf(3,...,n) = 7Ve(r''^..T'^"), Gr|:^^) (3; 4, n) = 0, Gr|:^li(3, n) = Tr (T'^^..r'^") <51^ , 

(36) 



Gri^f (3, + 1; J + 2, ...,n) = Tr(r°\..T''^-+i)(T'^^+^..r""):^, j = 3, - 2. 



As the next step, we describe the color decomposition of — it -I- -I- 17 amplitude [38|, where 
g is a massless quark. The color decompositions of the tree and the one-loop amplitudes read 



S'"'=(l,-,2t,3„4„5g 



1 



(37) 



Si-i°°P(lf,2i,3,-,4„5g) 5s 



A^c(r"^)^CBi + (T''=)^C^2 



^c(r"^),rc^3 + (T'^^).re54 



(38) 



In each of these one-loop color-ordered amplitudes we separate terms with the closed massless 
fermion loop, with the top quark loop and all other contributions 



i = 1,2,3,4. 



(39) 



The amplitudes and ^'5*^*'! can be written as linear combinations of primitive amplitudes. 
Those primitive amplitudes are illustrated in Fig. [T] Explicit expressions for color-ordered ampli- 
tudes bI'^^ and in terms of the primitive amplitudes are given in Ref. [11]; they can be 
found in the Appendix A. 
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Figure 1: Prototype parent diagrams for primitive amplitudes with four quarks, and a gluon for classes "a", "b", 
"c" and closed fermion loops (d). The gluon can be inserted in four possible ways into any of the prototype graphs, 
leading to four primitive amplitudes in each case. 



In Ref. [16j] , numerical results for spin-averaged squared amplitudes for the virtual corrections 
were given for a particular kinematic point. In Appendix B, we give numerical results for most of 
the helicity amplitudes for that kinematic point. We have checked the calculation in several ways. 
For example, our primitive amplitudes have correct 1/e poles 5lJ and are gauge invariant, once 
mass counter-terms are accounted for. We performed a diagrammatic check comparing amplitudes 
computed through unitarity and on-shell matrix elements with a similar, but independent, imple- 
mentation of the OPP procedure |35| which is applicable to individual Feynman diagrams. Finally, 
we have checked that when we compute the spin-averaged amplitude squared for the kinematic 
point considered in Ref. [l^, we find full agreement with their results. 

A well-known problem that occurs in the process of the reduction of one-loop diagrams to 
any basis of scalar integrals is the appearance of unphysical singularitiet|f| at the intermediate 
stages of the calculation. Those singularities cancel in the final result as a matter of principle 
but it has been proven to be difficult to achieve such cancellations in practice if reductions are 
performed numerically. Within the current calculation, we control potential numerical instabilities 
by requiring that the 1/e pole of a primitive amplitudes is calculated correctly to within 0(10"^). 
If this condition is not fulfilled, the program switches from double to quadruple precision and 
re-calculates the offending amplitude. The number of points that needs to be re-calculated is 
typically small, below a percent for partonic channels that we consider. 



3.3. Real emission processes and soft-collinear subtraction counter-terms 

Computation of the real emission corrections and subtraction counter-terms is an important 
part of any NLO computation. There are three generic processes that need to be considered: 
— >■ tt+gggg, — tt+qq+gg and — > tt+qq+q'q' . Once amplitudes for those processes are known, 
all partonic channels that contribute to the real emission correction to the hadronic production of 
tt + jet can be constructed. We employ color decomposition to exp ress those amplitudes through 
color-ordered ones and use Berends-Giele recurrence relations [52| to compute the latter. Below 
results for color decomposition of the relevant tree amplitudes are summarized. 

The color decomposition for the amplitude ^ tt + gggg has already been given in Eq. (j33p . 
The color decomposition for the amplitude ^ tt + qq + gg in double index notation for the gluon 



* Those unphysical singularities are caused by the so-called Gram determinants. 
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field reads [53 1 



Finally, the color decomposition for the six quark amplitude for non-identical quarks — ^ tt 
m + // reads 

A{t,t,q,qjj) = S^S;^qA(t,t,q,qJj)+S^_5^S^A(i,tJj,q,q) 

-W i^^^^^n + m^";^ - W^^Z^^ AitJ,q,qJ-,t) 



The amplitude for the identical quarks A{t, t, q, q, q, q) can be obtained from the previous 
amplitude by the anti-symmetrization procedure 

A{t, t, q, q, q, q) = A{t, t, q, q, f, f) - A{t, t, q, f, f, q). (41) 

We note that many of the amplitudes required for the computation of the real emission cor- 
rections are also employed for the computation of the virtual corrections within the unitarity 
framework. We use one and the same program that computes the on-shell tree amplitudes for 
complex or real external momenta, as needed for the virtual and the real corrections, respectively. 
We checked our calculation of the real emission matrix elements against Madgraph 5J| for all the 
partonic channels. 

As it is typical for NLO calculations, when we integrate the real emission matrix elements with 
tt and two additional partons in the final state over the phase-space constrained by the requirement 
that a tt pair and at least one jet is observed, we encounter infra-red and coUinear divergences. 
Those divergences are removed by the subtraction procedure. We use dipole formalism of Ref. [s^ 
extended to deal with the QCD radiation off massive particles in Ref. [H^. We note that in 
the original publications, the subtraction terms were integrated over full unresolved phase-space 
which is not very convenient. An optimization of the subtraction procedure, where subtractions 
are only performed if the event kinematics is close to singular, was suggested in Ref. [H^l- Dipoles 
integrated over the restricted phase-space including the case of massive emitters and spectators 
can be found in Ref. [23, 58, 5^ 6^. In the actual implementation of the subtraction procedure, 
we closely follow Ref. [58[ . We have checked that our results do not depend on the parameter that 
restrict the integration over the dipole phase-space; this is a useful way to check the consistency 
of the implementation of the subtraction terms. 



4. Results 

In this Section we present the results of the calculation of the NLO QCD corrections to tt pair 
production in association with one hard jet at the Tevatron and the LHC. We begin by comparing 
our results with that of Refs. 1^ 16 1. Then we present the results that include the top quark 
decays. 
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Figure 2: Top quark transverse momentum (a) and rapidity (b) distributions (MS) computed through NLO QCD 
for pp — >• ttj, at the LHC (y/s = 14 TeV) compared to the results of Refs. [H, (DUW). The renormalization 
and factorization scales are set to fi = rat- 



Input parameters and comparison with existing results 

For the comparison with Refs. [H, [l6|, we need to choose identical input parameters. To 
this end, we use the top quark mass mt — 174 GeV. Also, t he p arton distribution functions 
CTEQ6L1/6M have been updated since the results of Ref. 3, 3| appeared; therefore, for the 



comparison we employ their older versions. Finally, the fc^-clustering algorithm with Rij — 
\/(yi — yi)^~ (4'i ~ (t'i)'^ = 1 and i?x-weighted recombination scheme [6l| is used to define jets 
in Refs. [15l |l6[. The jet transverse momentum cut p±j — 50(20) GeV is employed for the 
LHC (Tevatron), respectively. When we use those parameters in the calculation, we obtain good 
agreement with Refs. [l^ Eal- We quote our results for the LHC {^/s — 14 TeV), with the 
factorization and the renormalization scales set to mt- For the total cross-section we obtain 

a^^'^ipp ^ ttj) = 375.8(1.0) pb, (42) 



which agrees well with the result 376.2 pb quoted in Refs. [15|, |16|. We note that the NLO 
QCD pp — ^ ttj cross-section at the LHC for the input parameters of Refs. [l^ [iBl was also 
reported in Ref. where the result 376.6 pb was found. We have checked that similar level 
of agreement - one percent or better - persists for cross-sections evaluated for other values of 
the renormalization and factorization scales, for both the Tevatron and the LHC. We have also 
verified that we reproduce kinematic distributions presented in Refs. As an illustration, 



our results for the transverse momentum and the rapidity distributions of the top quarks at the 
LHC, computed through NLO QCD, are compared with similar distributions from Refs. 1^ 3 
in Fig. [2] G ood agreement between the two results is apparent. We also confirm an observation 
m 

Refsnil,[il 

that the forward-backward asymmetry of top quarks in tt + jet production at the 
Tevatron is significantly reduced if NLO QCD effects are taken into account. For fi = mt and 
the jet transverse momentum cut of 20 GeV, we find the forward-backward asymmetry computed 



through NLO QCD to be -2.28%, in agreement with -2.27% reported in Refs. [15|, [lef. We 
discuss the forward-backward asymmetry in more detail below, when we report on our calculation 
that includes decays of top quarks. 

Having performed an extensive comparison with results of Refs. [isl. [l6l|. we turn to the discus- 
sion of observables that can be studied if decays of top quarks are accounted for in the calculation. 
We consider two primary scenarios - ttj as the signal process and as the background process to 
weak boson fusion production of the Higgs boson. It is important to consider these different sce- 
narios because we want to understand to what extent applied cuts affect the radiative corrections. 
We begin with the discussion of the pp{pp) — ttj production in the kinematic region which is usu- 
ally employed to study top quarks. We perform separate studies for the Tevatron and the LHC. 
We choose two center of mass energies -^/s = 7 TeV and y/s — lA TeV for the LHC and one center 
of mass energy -^/s = 1.96 TeV for the Tevatron. We use the following values for the top quark and 
the VF-boson masses mt = 172 GeV, mw = 80.419 GeV and employ CTEQ6L1 and CTEQ6.1M 
parton distribution functions in LO and NLO calculations, respectively. We employ the k± jet 
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Figure 3: Various distributions for the process pp — > (t — > I'^ub) + {t jlj2ti) +j at the Tevatron at leading (blue) 
and next-to-leading order (red) in perturbative QCD. The bands correspond to the choice of the renormalization 
and factorization scales fi = [mt /2, mt , 2mt] . We show distributions of the transverse momentum (a) and rapidity of 
the positron (b) , the total transverse energy H± (c) , the transverse momentum (d) and the rapidity (e) distributions 
of the fifth hardest jet and the missing transverse momentum (f). 

algorithm with i?y — 0.5 and the four-momentum recombination schem^. The couplings of the 
H^-boson to fermions are obtained from the Fermi constant Gp = 1.16639 • 10^^ GeV~^. We em- 
phasize that results reported in this paper are calculated with on-shell top quark decays at leading 
order in QCD; radiative corrections to the decays are not included. We take the CKM matrix to 
be the identity matrix. We consider the M^-boson in t — > Wb decay to be on the mass shell. With 
these approximations, we obtain the leading order top quark decay width Ft — 1.47 GeV. We set 
the width of the W boson to 2.14 GeV. With our input parameters, we find the following branch- 
ing ratios Br(VK+ l+v) = 10.6 %, Br(M^+ -> hadrons) = 64 %. We note that the hadronic 
branching fraction is slightly lower than the experimental value 67 %. 

4-.2. tt-|-jet production in lepton + jets channel at the Tevatron 

We begin our discussion with the Tevatron, ^/s = 1.96 TeV. The measurement of the pp — 
tt + jet at the Tevatron was performed by the CDF collaboration and there is a public note 
[Hi! that describes the details of the experimental setup and the results of the measurements. 



^ We note that Rij = 0.5 in the fc^ algorithm roughly corresponds to Rij = 0.4 in a typical cone algorithm [6^ . 
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For our purposes, we can take all the selection criteria for pp — >■ ttj from that reference. We 
consider events where top quark decays leptonically and the anti-top quark decays hadronically. 
We require that there are five or more jets in the event and that two of these jets are 6-jets. It is 
to be understood that a 6-jet is defined as a jet that originates from the tWb vertex. For pp — tij 
process, this leads to a misidentification of 5-jets, that originate from tt production, as non-6 jets. 
However, this is a minor issue since 6-jets in the production process appear infrequently. The 
lepton transverse momentum and the missing energy in the event should satisfy p^j > 20 GeV, 
^miss ^ 20 GeV. The jet transverse momenta are required to be larger than p±,j > 20 GeV and 
jets must be central \yj\ < 2. There is an additional cut on the transverse energy in the event 
H± — J2jP^-j +P±.e + E™^^^, needed to better discriminate against the background. It is required 
that Hj_ > 220 GeV. 

We first present the results for total cross-sections. We set the renornialization and factorization 
scales equal to each other and consider three values fi = [mt/2, m^, 2mt]. We obtain 

<^Lo{pp -^{t^ l+uh) + {t^jij2b)+j) = 36.9i?2-i fb; 

<JNLo{pP^ {t ^ l+vb) + {t^ji]2h)+j) = 33.61^:^ fb, (43) 

where the central value refers to /i = mt, the upper value to ^ — mt/2 and the lower value to 
/i = 2mt. The improvement in the scale stability at next-to- leading order is apparent. 

Our results for the kinematic distributions of the transverse momentum and the rapidity of the 
charged lepton the total transverse energy the transverse momentum and the rapidity of 
the fifth hardest jet and the missing transverse energy are shown in Fig.[3][a)-(f), respectively. The 
bands correspond to choices of the renormalization and factorization scales = [mt/2,mt,2mt\. 
We observe that the scale variation in rapidity distributions is in line with the total cross-section 
modification and there is little change in shape. However, for the transverse momenta distributions 
and the H±^ distribution, the situation is more subtle. There are clear differences in shapes between 
distributions computed at leading and next-to-leading order and, in addition, there are kinematic 
regions where the scale variation bands at LO and NLO do not overlajzH. In particular, this 
happens at the high end of the lepton transverse momentum distribution where, in addition, the 
scale variation at NLO actually exceeds the scale variation at leading order. Similar phenomenon 
occurs in the missing transverse energy distribution. On the contrary, the scale variation of the 
transverse momentum distribution of the fifth hardest jet is modest at high p_lj. 

It is interesting that some of the distributions exhibit peculiar shape distortions. We stress 
that these distortions depend in a very significant way on the choice of the renormalization and 
factorizations scales in leading order computations and that it is entirely possible that some of 
the changes that we observe in Fig. [3] can be accommodated by a kinematic-dependent choice 
of the renormalization scale. We will not pursue this topic further in the current paper; for a 



related discussion in the context of -I- jets production, see Refs. [6J, |65|, |66|, |67| . For our choices 
of scales, we observe that the charged lepton transverse momentum distribution and the missing 
transverse energy distribution become softer at NLO when compared to corresponding leading or- 
der distributions. We can understand this by observing that foru = m^, the transverse momentum 



distribution of the top quarks becomes softer at NLO as well [15|, |16[. On the other hand, fi = mt 
is almost the perfect choice of the renormalization scale for the transverse momentum distribution 
of the fifth hardest jet; this is probably just a numerical coincidence. While H± distribution does 
not show significant distortion at high H± , the (broad) peak of the H± distribution shifts to higher 
Hj^ values and there is a depletion at low H± values. 

Finally, it is clear from Fig. [3] that the positror0 rapidity distribution becomes much more 
symmetric at next-to-leading order. To quantify this, we compute the positron forward-backward 



®This statement is, of course, contingent upon the chosen range of scales. 

''We will refer to the positively charged lepton as the "positron" in what follows, but everything that is said 
applies to positively charged muons as well. 
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asymmetry 

For the set up described above, we find the positron asymmetry to be strongly reduced at next- 
to- leading order, similar to the top quark asymmetry 1^ T^. We obtain 



A'^V - -5.05 %, = -0.5 %. (45) 

At first sight, the strong reduction in the positron forward-backward asymmetry observed in 
pp — >■ ttj process at NLO is rather puzzling and worrisome since it suggests a breakdown of the 
perturbative expansion for this quantity. We will now argue that i) these worries are unfounded; 
ii) the large reduction in the forward-backward asymmetry in pp — >■ ttj at NLO is natural; and 
in) the NLO asymmetry is, most likely, stable against yet higher order corrections. To build up 
the argument, we consider the limit of the low jet transverse momentum cut p±j — )■ 0. 

We notice that, by requiring additional jet in the final state and by taking p±j to be small, 
we introduce two types of degrees of freedom: soft degrees of freedom, controlled by the jet 
transverse momentum cut p±j and hard degrees of freedom, controlled by the mass of the top 
quark. The crux of our argument is that the asymmetry in tij can be generated by both soft and 
hard degrees of freedom but those mechanisms appear at two consecutive orders in perturbation 
theory. The "soft" asymmetry appears at leading order while the "hard" asymmetry appears 
at next-to- leading order. As evident from their very different dependence on p±j, these two 
mechanisms are unrelated. Hence, the dramatic change in the asymmetry at next-to-leading order 
should not be taken as an indicator that the perturbative expansion for this observable breaks 
down. 

We now explain this argument in detail. We discuss the top quark asymmetry rather than 
the positron asymmetry, for simplicity. In the limit of small p±j, the leading order asymmetry 
is generated when the soft gluon between initial (qq) and final (ti) states is exchanged. Since 
the interference diagrams are non-singular in the collinear limit, the difference of the forward and 
backward cross-sections at leading order can only be proportional to a single logarithm of the 
infra-red cut-off 

la{yt > 0) - <j{yt < 0)]lo « <JA In ( ^) , 2/t - ^ In (^^^) ■ (46) 

\P±,j/ 2 \Et-pt,zJ 

Here, a a is some quantity with the dimension of the cross-section and yt is the top quark rapidity. 
On the other hand, the total cross-section that appears in the denominator in the definition of the 
asymmetry Ea. (|44p depends on the double logarithm of the jet transverse momentum cut 

a{yt > 0) + a{y, < 0) ^ ^ In^ ( ^\ a,,. (47) 



In Ea. (|Tf|) . (Ttf is the production cross-section for pp — > ti. As the result 



iLO/„ ^ _ (^{yt > 0) - a{yt < 0) _ ( rut 



^^^yf^^i! a(yt > 0) + a(yt < 0) \p^^J ^ ' 

Hence, we see that at leading order in perturbative QCD, the asymmetry in pp — > ttj is gener- 
ated by soft degrees of freedom and that this asymmetry decreases with the decrease in the jet 
transverse momentum cut 

At next-to-leading order, the asymmetry can still be generated by soft exchanges. In such a 
case, it is natural to expect that moderate correction to the leading order asymmetry Eg. pS)) is 
generated. However, it is interesting that at next-to-leading order, a new mechanism for generating 
the asymmetry in pp — > ttj appears. Indeed, we can use hard degrees of freedom to generate the 
asymmetry, in much the same way as it is generated in the inclusive pp — tt process. In addition, 
we can use the emission of an additional soft and collinear jet from the initial state to provide 
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Figure 4: Various distributions for the process pp {t l^vb) + (f — > l~vb) + j at the LHC {^/s = 7 TeV) 
at leading (blue) and next-to-leading order (red) in perturbative QCD. The bands correspond to the choice of the 
renormalization and factorization scales fi = [mt/2,mt,2mt]. We show distributions of the transverse momentum 
(a) and rapidity (b) of the positron, the total transverse energy H± (c), the transverse momentum distribution of 
the third hardest jet (d) and the invariant mass of the two leptons (e). 



conventional double logarithmic enhancement of this "hard" asymmetry. This mechanism leads to 
the following forward-backward cross-section difference in Eq. ([^^ as 

aiy, > 0) - a{y, < 0) ^ In^ ( ^] A,, a,,, (49) 

where A^^ is the pp ti inclusive asymmetry. Taking the ratio of Eg . (1491) and Eq.(j47l), we find 

A^ij" « Au. (50) 

The full forward-backward asymmetry at next-to-leading order is given by the sum of the two 
mechanisms 

A^o«A?g*(p^,,)+A^'^. (51) 
Following the preceding discussion, we estimate 

A?f(p^,,)^ALO(p^.,)^ln-i(^^), (52) 
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Numerically, the "soft" asymmetry at p^j = 30 GeV is A]'^{30 GeV) -7% and the "hard" 
asymmetry is A^^ ~ 5% (Hsj , so that the significant reduction in the leading order pp — >■ ttj 
asymmetry is observed, once the NLO QCD corrections are accounted for. A peculiar consequence 
of this argument is that, in the limit of a very small jet transverse momentum cut, the top quark 
forward-backward asymmetries in the pp ttj process and in the inclusive pp — >■ ti process 
coincide 

lim Atij {p±,j) Ati, (53) 

because the "soft" component of the asymmetry vanishes as the inverse logarithm of the cut on 
the jet transverse momentum. 

We stress that the observation by Dittmaier, Uwer and Weinzierl of the large NLO QCD 
corrections to the top quark forward-backward asymmetry in pp — > ttj [illil is important outside 



the context of that computation. Indeed, since the asymmetry in pp — tt appears at NLO and since 
NLO is the highest order in the perturbative expansion available for that process, corrections to 
the inclusive asymmetry are not known. On the other hand, because the asymmetry in pp — ;> ttj 
appears already at leading order, a NLO QCD computation for pp — ^ ttj gives corrections to 
the leading order forward-backward asy mmetry. The fact that these corrections turn out to be 
nearly 100%, as discovered in Refs. la llfll: leads to doubts about the robustness of the existing 



predictions for the inclusive asymmetry in the pp — tt process. Of course, the issue of robustness 
is rather iniportant given the existing discrepancy between the measurement of the asymmetry in 
pp — d, Q and the theoretical prediction [68] . 

We believe that our argument supports the robustness of the theoretical prediction of the in- 
clusive asymmetry Aft. In fact, as we have explained, there are two physically distinct mechanisms 
that generate the asymmetry in pp — ttj. These mechanisms are related to the existence of two 
types of degrees of freedom, which can be clearly separated by lowering the cut on the jet trans- 
verse momentum. The only peculiarity about the asymmetry is that the two mechanisms for the 
asymmetry generation do not appear at the same order of the perturbative expansion which leads 
to apparent problems with the interpretation of its convergence. Needless to say that we do not 
see any other mechanism that can start contributing to the asymmetry at next-to-next-to-leading 
order (NNLO) and beyond. Therefore, we believe, that the NLO prediction for the asymmetry in 
pp — ttj is robust. Similarly, since inclusive pp — tt production process is only sensitive to hard 
degrees of freedom, no new mechanism for the asymmetry generation can appear at NNLO and 
beyond. Therefore, we believe that the existing prediction for the top quark forward-backward 
asymmetry in pp — )■ ti is robust as well. 

4.3. tt-|-jet production in dilepton channel at the LHC 

We now turn our attention to the LHC, y/s — 7 TeV. We consider dilepton channel and require 
three or more jets with p±,j > 50 GeV, p±,b-jct > 20 GeV. Our definition of the 5-jet is explained 
above. We also apply the following cuts on the lepton transverse momentum p±^i > 20 GeV and 
the missing transverse momentum p_L,miss > 40 GeV. We begin by showing the cross-sections for 
the three values of the factorization and renormalization scales /i = [mt/2,mt,2mt\. We find 

cTLoipP {t^ l+iyb) +{t^ l-9h)+j) = 229.9^7^^/ fb; 

CTNLobp^ it l^J'b) + [t^ l-Dh)+j) = 256.51^^:^ fb, (54) 

where the central value refers to ^ — nit, the upper value to /i = mt/2 and the lower value to 
/i = 2mt. Similar to the Tevatron case, the NLO cross-section is very stable against changes in 
the renormalization and factorization scales, in contrast to the leading order result. 

Our results for kinematic distributions are displayed in Fig. 21 where we present the transverse 
momentum and the rapidity distributions of the positively charged lepton, the distribution of the 
total transverse energy H^, the transverse momentum distribution of the third hardest jet and 
the distribution of the invariant mass of the two leptons. We note that for some choices of the 
renormalization and factorization scales, the NLO results for the lepton transverse momentum 
distribution become negative at high p±{l^)^ 400 GeV. This indicates that /i = mt/2 is too low a 
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Figure 5: Kinematic distributions for the process pp —> {t —> l^vb) + (t — > l~vb) + j at the LHC (v^ = 14 TeV) 
with weak boson fusion cuts. The renormalization and factorization scales are set to = mt- Distributions at 
leading (blue) and next-to-leading (black) order in perturbative QCD are shown. We show distributions of the 
transverse momentum (a) and rapidity of the electron (b), the invariant mass of the electron and the positron (c), 
the transverse momentum of the leading jet (d). 



scale for this observable and that there is large residual scale dependence in this observable even at 
next-to-leading order. The 77j^-distribution is shifted to higher values, leading to positive correc- 
tions in the high-_ff_L tail. The transverse momentum distribution of the third jet becomes softer 
and the lepton invariant mass distribution is only marginally affected. Many of these features can 
probably be understood by performing a leading order computation with the kinematic-dependent 
choices of the scales. As we already mentioned in the context of the Tevatron discussion, we do 
not pursue this topic in the present paper and reserve it for future work. 

Our last example concerns NLO QCD corrections to pp — > ttj process in the kinematic region 
relevant for the Higgs boson searches in weak boson fusion pp — > jiH jjW^W^ production 
channel, with leptonic decays of the IV-bosons. It was shown in Rcf. [69| that the pp — > ttj process 
is the largest background to the weak boson fusion signal. For our calculation, we take the LHC 
energy to be =14 TeV. The event selection is based on the following cuts '69|. We define jets 
with the usual k±_ algorithm, Rij = 0.5 and the transverse momentum cut p±,j > 20 GeV. The 
two jets with the highest transverse momentum (tag jets) are required to satisfy p^"* > 40 GeV 

(2) 

and p\ > 20 GeV. Those jets should be in the opposite hemispheres yij/2 < and widely 
separated in rapidity \yi — 2/2I > 3.0. The invariant mass of the two tag jets should be large 
rrij-^j.^ > 550 GeV. The two leptons are required to have large transverse momenta p±,i > 20 GeV 
and central rapidities \yi\ < 2. In principle, there are other cuts that are imposed to separate the 
weak boson fusion signal but we do not employ those cuts in what follows. For the cross-sections 
at 11 — mt we find 

<JLo{pp -> (t ^ l^i^b) + (t ^ l~Db) + j) = 139.6 fb; 

(7NLo{pp -^{t^ l+vb) +{t^ rDh) + i) = 177.9 fb. (55) 

The next-to-leading order cross-section exceeds the leading order cross-section by twenty-five per- 
cent which is not unusual. In Fig. [S] we show a number of kinematic distributions for the weak 
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Figure 6: Kinematic distributions for the process pp — > — > l^ub) + (f — > l~Pb) + j at the LHC = 14 TeV) 
with weak boson fusion cuts. The renormalization and factorization scales are set to fi = mt- Distributions at 
leading (blue) and next-to-leading (black) order in perturbative QCD are shown. We show distributions relevant 
for the discrimination between the signal and the background. In particular, we display distributions of the relative 
azimuthal angle of the electron and the positron (a), the my variable (b) (see text), the rapidity of the leading jet 
(c), the rapidity of the second leading jet (d) and the relative rapidity of the "veto" (third hardest) jet (e). For the 
definition of the relative rapidity of the veto jet, see text. 



boson fusion cuts. In particular, lepton transverse momentum and rapidity distributions, distri- 
bution of the dilepton invariant mass and the transverse momentum distribution of the hardest 
jet are displayed. In Fig. |6l we show kinematic distributions that can be used to discriminate 
between the tt + jet background and the weak boson fusion signal. We display the distribution 
in the relative azimuthal angle of the two leptons, the distribution of the approximate transverse 
mass of the two W^-bosons tot and some rapidity distributions. The definition of the TOy-variable 
is given in Ref. ^] 



rriT = 




In this formula, p_L,miss is the missing transverse momentum which we associate with the vector 
sum of the momenta of the two neutrinos and Ej_^^iss — \J p\ miss + ''^f+i- ■ explained in 
Ref. [i^ , rriT is a good approximation to the actual transverse mass of the Higgs boson, provided 
that the Higgs mass is betwen 160 and 200 GeV. A glance at Figs. [5l [6] suggests that shapes of 
lepton p± and rapidity distributions and the dilepton invariant mass distribution are not strongly 
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affected by the radiative corrections for fi — mt- Tlie distribution of the transverse momentum of 
the leading jet becomes harder. The shape of niT distribution is not changed by the NLO QCD 
corrections but there appears to be a shape change in the distribution of the relative azimuthal 
angle of the two leptons. As the result, the A(j> distribution at next-to-leading order is less peaked 
at A(f) sa TT, compared to the leading order result. We also plot two rapidity distributions. The 
rapidity distribution of the hardest jet at NLO (Fig. EJc)) is well reproduced by the re-scaling of 
the leading order distribution. On the other hand, the NLO rapidity distribution of the second 
hardest jet, shown in Fig.|6l(rf), is shifted to larger absolute values and the shape of the distribution 
changes. The distribution of the relative rapidity ?/roi = 2/veto — 0.5 (jjj-^ + yj^) of the "veto jet" is 
also shown. Fig. Elje). The veto jet is defined as the jet between the two tagged jets. The shape 
of j/ioi distribution does not change from leading to next-to- leading order. 



5. Conclusions 

We describe the calculation of the NLO QCD corrections to the production of a tt pair in 
association with a hard jet at the Tevatron and the LHC. The one-loop virtual corrections are 
computed using the framework of generalized Z?-dimensional unitarity. We present numerical 
results for one-loop helicity amplitudes for the processes — tiggg and ttqqg. 

Availability of helicity amplitudes enables us to include decays of top quarks, in the narrow 
width approximation, at almost no additional cost. We do that in leading order in perturbative 
QCD. We account for all the spin correlations between top quark decay products exactly and 
include constraints on top quark decay products in the computation. We emphasize that these 
features may be important for the detailed comparison between the theoretical predictions for 
pp(j)p) tt + jet and results of experimental measurements as well as for the the understanding 
of the tt + jet process in the kinematics relevant for New Physics searches. 

We consider production of the tt pair in association with one hard jet both at the Tevatron 
and the LHC, for realistic cuts on the final state leptons, jets and missing energy. For the LHC, 
we also investigate NLO QCD corrections to tt + jet process in the kinematic region selected by 
the weak boson fusion cuts 69]. For all the cases considered, we find reduced dependence on the 
renormalization and factorization scales. This is true for total cross-sections and main regions 
of kinematic distributions. For tails of kinematic distributions, the situation is more subtle since 
next-to-leading order and leading order scale-dependence bands are often similar and may not 
even overlap. 

We confirm findings of Refs. 1^ 1^ that the top quark forward-backward asymmetry in the 



pp — >■ ttj process receives very large corrections at next-to-leading order. We explain the origin 
of these large corrections and argue that they should have been expected. Since there is no 
mechanism, that would lead to large corrections to top quark forward-backward asymmetry in 
pp — > ttj process in yet higher orders, we consider the NLO QCD prediction for the asymmetry 
as quite reliable. We note that large NLO QCD corrections to the asymmetry in pp — > ttj process 
are particular to the presence of a jet in the final state; the enhancement mechanism does not 
work for the inclusive asymmetry in pp — )■ tt. Therefore, we do not expect that missing higher 
order corrections are the reason for the significant discrepancy between the QCD prediction for 
the forward-backward asymmetry (ssj and the experimental results 0, 0] . 

The calculation reported in this paper opens up a possibility to study production of a tt pair in 
association with a hard jet at hadron colliders in a realistic manner. Further improvements may 
require the inclusion of the QCD corrections to top quark decays. In addition, there is significant 
interest in other processes where a ti pair is produced in an association with either a vector boson 
or the Higgs boson l^, 2^ . It will be definitely interesting to refine NLO QCD predictions for 



those processes to account for the decays of top quarks. 
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Appendix A. 



In this Appendix, we give explicit expressions 38| for color-ordered amplitudes Bf'^ and 



^ji/2,9(t)] terms of the primitive amplitudes. The different classes of the amplitudes are shown 
in Fig. [TJ For the amplitudes with closed fermion loop we find 

There are three classes of primitive amplitudes that we need to consider for four-quark processes 
without closed fermion loops 

^[1] ^^[l],a^^[l],6^^[l],c^ (A.2) 

they are shown in Fig. [Ija)-(c). Amplitudes from each class are written as linear combinations of 
primitives amplitudes. 

For the class "a" , we find 

M'''" - (^1-^)aW'"(1,,2„3„4„5,)-^(-41''^(1,,5„2„3„4,) (A.3) 

-A™'" (If, Sg, 2i, 4„ 3,-) - A™'" (If, 2t,5g, 3,-, 4,) - A^l^'"" (If, 2^,63,4,, 3,-) 
-A™'" (If, 2f, 3,-, 5g, 4,) + A™'" (If, 2f , 4„ 5g, 3,-) + A™'" (If, 2* , 4„ 3,-, 5^)) , 
^[i],a ^ ^li],a 53, 4,, 3,-) + A™^" (If, 2t, 4g, 5^, 3,-) + A™''^ (If, 2f, 4„ 3,-, 5g) 

(4"'" 5s, 2f, 3,, 4,) + aW'"^ (If, 5g, 2f, 4„ 3,)) , (A.4) 

4'''"^ = (^l-^)41^"(lf,2f,5„3„4,)-i^(-41''^(lf,5„2f,3„4,) (A.5) 

-A™'" (If, 53, 2f, 4„ 3,-) - A™'" (If, 2t, 3,-, 53, 4,) + A™'" (If, 2^,4,, 5,,, 3,-) 
-4^1'" (If, 2t, 3,-, 4„ 5g) - A™^'^ (If, 2f , 4„ 3,-, 63) + A™'" (If, 2*, 5g, 4„ 3,-)) , 
4"'" = -4"'" Sg, 2f, 4„ 3,-) - A™^" (If, 2^,63, 4„ 3,-) - A™'" (If, 2^,4,, 3,-, 5g) 

-772 (4"'" 2t, 3,, 5„ 4,) - aW'"^ (If, 2t, 4„ 5„ 3,)) . (A6) 

For classes "b" and "c" we obtain 



M"'' = ^4''''(lt:5„4„3„2f), 



4"'" = (4''"' (It- 5g, 4„ 3,-, 20 -f 4^1^' (If, 4„ 5g,3g, 2*) + A™'" (If, 4„ 3,-, 5g,2t) ) , 



4'''" = ^4''''(lt:4„3„5„2f), 



M'^'" = -4"'' (1*: 5s, 4„ 3,-, 2,) - A™'" (If, 4„ 3g, 5^, 2f) - A™'" (If, 4„ 3,-, 2*, 5^) 

l-^)4''''(lt,4„5s,3,,2f), (A7) 



M'''^ = ^4'''(lt,5s,4„3„2f), 

4"'' - -4''' (lt> 59-4,, 3,-, 2f) - 4''= (If, 4„ 53, 3,-, 2f) - 41'^ (If, 4„ 3,-, 5g, 2*) 
1 



l-]^)4'''(lt,4„3„2t,5g). 
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l2 4'''(lt,4„3„5g,20, (Ai 
(4"'' (It, 5g, 4„ 3,-, 2t) + 4''" (It, 4„ 3,-, 5g, 2*) + 4^''= (1,-, 4„ 3,-, 2^,5^)) . 
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Appendix B. 

In the Appendix we give sample results for the helicity amplitudes for the kinematic point 
considered in Ref. [11] . The momenta for the reaction ab ^ ti + c read (in units of GeV) 

Pa = (500,0,0,500), 
Pb = (500,0,0,-500), 

Pt = (458.53317553852783, 207.0255169909440, 0,370.2932732896167), 

Pt = (206.6000026080000, -10.65693677252589, 42.52372780926147, -102.39982104210421085), 
Pc = (334.8668220067217, -196.3685802184181, -42.52372780926147, -267.8934522475083). 

We present results for the ratio of one-loop helicity left-primitive amplitudes to the corresponding 
tree primitive amplitudes in Tables B1-B4. Those results do not include external wave function 
renormalization constants and the coupling constant renormalization. However, mass counterterms 
are included to obtain a gauge- invariant result. We define 



rt^l(zi,i2,i3,i4,i5) 



1 A^l\ii,i2,i3,i4,i5) 
cr A*™<=(ii,i2,«3,i4,i5)' 



(B.l) 



The one-loop scattering amplitudes presented below are computed for definite helicity states 
of external quarks and gluons. For a gluon with momentum p = i?(l, sin0cos0, sinsin^, cos0), 
we define polarization vectors as 



^t{p) — ~~j= (0, cos cos (/i =F j sin0, cos sin ± i cos 0, — sin i 
v2 



(B.2) 



We use the Dirac representation of the 7- matrices. For a fermion with the on-shell momentum 
y = {E,px,Py,Pz), P^ — 'm? , and helicity A = ±, we use the spinors 



u+{p) = 







\ {pa; + ipy) / VE~frn J 

( PzHE + m \ 
(Px + ipy)I^E + m 



V 



yjETm 




u-{p) 



v-{p) = 







y/E- 



[Px - ipy) ly/E + m 
V -PzHE + m / 

( {Px - ipy) /y/E + m ^ 
—pz/y/E + m 




V 



(B.3) 



(B.4) 



\/E + m I 

The case of massless leptons is obtained by taking m — t- limit in the above formulas. 



23 



Helicity amplitude 






e 


2+,3-,4-,5+) 
rW(l+,2+,3-,4-,5+) 


-3.000000 


6.86872922 -6.28318531i 


-0.11654501 + 0.743827431 
8.70123544 + 3.769928171 


3-,2+,4-,5+) 
rW(l±,3-,2+,4-,5+) 


-2.000000 


5.61128078 -3.14159266i 


0.13444688 - 0.829874441 

7.23821384- 1.355966621 


>l'-^(ir,3-,4;,2+ 5+) 
rW(l±,3-,4-,2+ 5+) 


-1.000000 


4.47200426 - 6.28318345i 


-0.01014177 + 0.076362661 
6.07761708 + 8.884515241 


A*-^(lt,3-,4-,5+ 2+) 
rW (1^3-, 4-, 5+ 2+) 


0.000000 


2.73111593-3.223888471 


-0.00776010 + 0.009684361 
62.01806062+ 113.318111241 



Table B.l: Examples of A-jj amplitudes contributing to the process 9334 — > ttg^. The tree level amplitudes A*' 
axe given in units of (100 GeV)~^. 



Hcllclty amplitude 


1/e^ 


1/e 


6" 


A*>--(lt,2+.3,-,4-,5+) 
r[V2w](i+ 2+ 3-,4-,5+) 


0.000000 


0.00000000 


-0.11654501 + 0.743827431 
-0.06796975- 0.261811301 


A ''°°(lj , 5g , 2( , 3g , 4g ) 

rlV2..](i+ 5+ 2+ 3-,4-) 


0.000000 


0.00000000 


0.00077387-0.123091591 

-0.11916673 - 0.056442551 


A*>--(1+ 2r,3-,4-,5+) 
r[VV](i+ 2+ 3-,4-,5+) 


0.000000 


0.00000000 


-0.11654501 + 0.743827431 
-0.10289728 - 0.068883691 


A ''®*'(lj , 5g , 2( , 3g , 4g ) 

r[V2,*l(i+ 5+ 2+ 3-,4-) 


0.000000 


0.00000000 


0.00077387- 0.123091591 
-0.13739705 + 0.064723321 



Table B.2; Examples of A^^'^' amplitudes contributing to the process 3394 — > ttg5- First and second rows correspond 
to closed fermion loops with masslcss fcrmions; third and fourth rows correspond to closed fermion loops with 
massive top quarks. The tree level amplitudes J^""^" are given in units of (100 GeV)~^. 
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Helicity amplitude 


1 /^a 




^0 
e 


/I tree 


1 + 0+ /I + 


'5g) 






—0.03342616— 0.15320512i 


[11 a 




'5g) 


— z.UUUUUU 




D.ooioDozz + U.ooUoooUoi 


A tree 


1 + c— 9+ 


4 
I ^9 ; 






0.01451328 — 0.078627491 






. 4^ 1 


— i.UUUUUU 


4.y /zozyi / 


O.DioozOlz — i.DD4yyoiZl 


A tree 
A 




4 






—0.00409831 — 0. 014940841 


rii a 


.If" ) 2^'', 5g , 


4 


— Z.UUUUUU 


y.oooiioiy — o.i4ioyzDDi 


— z.UdoDDdoz + o.U/UzyoUol 


A tree 
/i 


'1+ 9+ 


>4g ) 






0.02301119 + 0. 246773451 


f ^ ' 


'1+ 9+ o- c- 




1 nnnnnn 




fi fi9nnnnn'i — 9 "^491 797'^i 

U.UZUUUUUO — Z.O'±Z± 1 Z 1 Ol 


^tree 


'1+ 2+ f)~ 4+ 








n ni SQ1 9SS n 9'?i s'?9fii i 
— u.u±oy±zoo — u.zo±oozD±i 




'1+ 9+ K— A + 




—2 flonnno 


8 747Q2056 — 3 1415Q267i 


9982.3846 + 2 119874471 


^tree 


'1+ 2+ 4+ 3~ 


.Og ; 






nnin4i4Q7 n no'^f^fiS'^'^^ 




'1+ 9+ A+ Q~ 
.J-t )«3g 




—2 nnnnnn 


Q Q75172m — 3 1415Q2fi7i 


—2 81 948074 + 3 4865n797i 


^tree 


2^ 4+ 5~ 

k-^t '"'t '*9 '^9 


'"*9 ) 






U.UzoUiiiy -r U.Z4D/ / O401 




'1+ 9+ A+ K — 

,t)g 


.<5g j 


—1 nnnnnn 


8 2QQ3n754 




^tree 


k t '"g ' ■^t > ^9 


; ^g ; 










'l+,5~,2;^ 4+ 


3^) 


-i.onnnnn 


8.29930751 


-4.93230098- 1.649745681 


^tree 


'1"*" 5~ 4+ 3~ 

k^t '^9 '^9 '"^9 


2+) 








/ L J 1 


1+ ^~ A+ 

'°9 '^9 '-^9 


9+^ 


n nnnnnn 


2 731 1 1 fil 2 — 3 2238884 7i 


8 8731 8696 + 3 n5n82436i 


^tree 


'^ + 4+ 5~ 3~ 

.J-i lOg j-Jg 


2+^ 








/ I J ' 1 


1+ 4+ ^~ 

'^q '^9 '"^9 


9+^ 


n nnnnnn 


2 731 1 1 5Qfi — 3 22388845i 


8 4846321 8+2 477nn393i 


^tree 


'1^ 4+ 3^ 5^ 


2+) 










1+ /!+ Q- 

if ,4, ,-J9 ,0, 


9+^ 


n nnnnnn 


2 73111596 — 3 22.388844i 


6 86621945 — 1 314n79121 


^tree 


'1+ 4+ 0+ 


^5g) 






— U.Ui40iozo -r U.U r oOZ r 4yi 




1+ 4+ 3" ■)+ 




-l.OOOOOO 


4.472()r71() - ().28323(i23i 


1().()()12")")")7 + <S.()31()81871 


^tree 


f-l+ 9+ o- c- 


4+) 






0.02301119 + 0.246773451 




1+ 9+ o- c- 
^f ' t ' "-^9 ' ^q 


4+) 


-2.onnnnn 


2.28510399 


2.72655473 + 4.843278711 


^tree 


,2j ,5g ,3g 


4+) 






-0.00409831 - 0.014940841 




1 + 9+ c- o- 
f ' t ' 9 ' <? 


4+) 


-I.onnnnn 


1.99739996 -3.14159265i 


4.68165727+ 2.921843371 


^tree 


'1+ 9+ '^^ 
k f ' g ' ■ <? 


4+) 

■ ^9 7 






0.01451328- 0.078627491 




1( ! 5^ , 2j , 3^ 


4+) 


-1.000000 


1.99739996- 3.141592661 


0.56499892 + 6.275017031 


^tree 


'1 + 9+ 4+ 


I ) 






-0.03342616 - 0.153205121 




1+ 2+ 3~ 4+ "i""! 

' ' "^9 ' ^q ' ^9 / 


-1.000000 


1.99739996-3.141592651 


2.76146380 + 3.498094871 



Table B.3: Examples of A^' amplitudes contributing to the process gs + 94 -> tt + ps. The tree level amplitudes 
A*''«e are given in units of (lOOGeV)-i. 
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i/fc 


1/fc 




[i/2,q]^+ 2+ 3r ,4+,5-) 


0.000000 


-0.66666667 


-0.03342616 - 0.15320512i 
0.14254796 - 2.21798474i 


^*-<=(lt 5;,2+,3r,4+) 
^[1/2,91(1+ 5-2+, 3-, 4+) 


0.000000 


-0.66666667 


0.01451328- 0.078627491 
1.22048886 - 2.094395101 


^*■■-(l^2r.5-,3^.4+) 

T)- 4+) 


().()()()()()() 


-().()()()()()()()7 


-0.00409831 - 0.014940841 
0.41410438 - 2.;-)22379711 


^tiooj-2_ ^ 2j , 3f^ , 5g , 4+) 
[i/2,q]/i+ 2+ 3- 5- 4+) 


0.000000 


-0.66666667 


0.02301119 + 0.246773451 
0.48191832 - 2.094395101 


r[i/V](l+,2+ 3r,4+ 5-) 


0.000000 


-0.66666667 


-0.03342616- 0.153205121 
-0.48762813 - 2.075145021 


, 5g , 2j , 3q , 4+) 
r[i/2.t](i+,5-,2+ 3r,4+) 


0.000000 


-0.66666667 


0.01451328 - 0.078627491 
1.08428102 - 2.082376841 


, 2( , 5g , 3^ , 4+) 
r[V2.*](i+ 2+ 5-,3r,4+) 


0.000000 


-0.66666667 


-0.00409831 - 0.014940841 

-0.26372154 - 2.569864491 


^''''^'^(1^ , 2( , 3g , 5g , 4+) 
r[V2.*l(i+,2+ 3r,5-,4+) 


0.000000 


-0.66666667 


0.02301119 + 0.246773451 
0.00939493- 1.972326581 



Table B.4: Examples of A^''^' amplitudes contributing to the process 93+94 -> tt + 35- The tree level amplitudes 
are given in units of (100 GeV)"!. 



26 



References 

[1] F. Abe et al. [CDF collaboration], Phys. Rev. Lett. 74, 2626 (1995). 
[2] S. Abachi et al. [DO collaboration], Phys. Rev. Lett. 74, 2632 (1995). 

[3] The Tevatron Electroweak fitting group, for the CDF and DO collaborations. larXiv:0803. 16831 
[hep-ex] . 

[4] A. Beretvas et al. [CDF collaboration] . larXiv:0710.4983[ hep-ex] . 

[5] V. M. Abazov et al. [DO coUaboration] , Phys. Rev. Lett. 98, 041801 (2007). 

[6] A. Abulencia et al. [CDF collaboration], Phys. Rev. D 75, 052001 (2007); Phys. Rev. Lett. 
98, 072001 (2007); Phys. Rev. D 75, 031102 (2007). 

[7] V.M. Abazov et al. [DO collaboration]. larXiv:071 1.00321 [hep-ex]. 

[8] T. Aaltonen et al, Phys. Rev. Lett. 101, 202001 (2008). 

[9] V. M. Abazov, et al, Phys. Rev. Lett. 100, 142002 (2008). 

[10] M. Mangano Eur. Phys. J. C 59, 373 (2009). 

[11] P. Nason, S. Dawson and R. K. Elhs, Nucl. Phys. B 327, 49 (1989) [ Erratum -ibid. B 335, 
260 (1990)]. 

[12] W. Beenakker, W.L. van Neerven, R. Menge, G.A. Schuler and J. Smith, Nucl. Phys. B 351, 
507 (1991). 

[13] M. L. Mangano, P. Nason and G. Ridolfi, Nucl. Phys. B 373, 295 (1992). 

[14] S. Frixione, M. L. Mangano, P. Nason and G. Ridolfi, Phys. Lett. B 351, 555 (1995). 

[15] S. Dittmaier, P. Uwer and S. Weinzierl, Phys. Rev. Lett. 98, 262002 (2007). 

[16] S. Dittmaier, P. Uwer and S. Weinzierl, Eur. Phys. J. C 59, 625 (2009). 

[17] G. Bevilacqua, M. Czakon, C. G. Papadopoulos and M. Worek, [arXiv: 1002.4009 [hep-ph]. 

[18] Peng-Fei Duan, Wen-Can Ma, Ren- You Zhang, Liang Han, Lei Guo, Shao-Ming Wang, 
Phys. Rev. D 80, 014022 (2009). 

[19] A. Lazopoulos, T. McElmurry, K. Melnikov and F. Petriello, Phys. Lett. B 666, 62 (2008). 

[20] W. Beenakker et al, Phys. Rev. Lett. 87, 201805 (2001); W. Beenakker et al Nucl. Phys. 
B 653, 151 (2003); S. Dawson, L. H. Orr, L. Reina and D. Wackeroth, Phys. Rev. D 67, 
071503 (2003); S. Dawson, C. Jackson, L.H. Orr, L. Reina, D. Wackeroth, Phys. Rev. D 68, 
034022 (2003). 

[21] A. Bredenstein, A. Denner, S. Dittmaier and S. Pozzorini, Phys. Rev. Lett. 103 (2009) 012002 
[arXiv:0905.0110 [hep-ph]]. 

[22] A. Bredenstein, A. Denner, S. Dittmaier and S. Pozzorini, JHEP 1003, 021 (2010) 
[arXi v: 1001. 40061 [hep-ph]] . 

[23] G. Bevilacqua, M. Czakon, C. G. Papadopoulos, R. Pittau and M. Worek, JHEP 0909, 109 
(2009). 

[24] W. Bernreuther, A. Brandenburg, Z. G. Si and P. Uwer, Phys. Lett. B 509, 53 (2001). 

[25] W. Bernreuther, A. Brandenburg, Z. G. Si and P. Uwer, Phys. Rev. Lett. 87, 242002 (2001). 



27 



[26] W. Bernreuther, A. Brandenburg, Z. G. Si and P. Uwer, Int. J. Mod. Phys. A 18, 1357 (2003). 

[27] A. Brandenburg, Z. G. Si and P. Uwer, Phys. Lett. B 539, 235 (2002). 

[28] W. Bernreuther, A. Brandenburg, Z. G. Si and P. Uwer, Nucl. Phys. B 690, 81 (2004). 



[29] W. Bernreuther, A. Brandenburg, Z. G. Si and P. Uwer, arXiv:hep-ph/0410197 
[30] K. Melnikov and M. Schulze, JHEP 0908:049 (2009). 



[31] W. Bernreuther and Z.G. Si. 'ar Xiv: 1003.3926 [hep-ph]. 

[32] W. T. Giele, Z. Kunszt and K. Melnikov, JHEP 0804, 049 (2008) [arXiv:0801 .22371 [hep-ph]; 

[33] R. K. Elhs, W. T. Giele, Z. Kunszt, K. Melnikov, Nucl. Phys. B 822, 270 (2009). 

[34] R. K. Ellis, W. T. Giele and Z. Kunszt, JHEP 0803, 003 (2008). 

[35] G. Ossola, C. G. Papadopoulos and R. Pittau, Nucl. Phys. B 763, 147 (2007). 

[36] Z. Bern, L. Dixon and D. Kosower, Nucl. Phys. B 437, 259 (1995). 



[37] W. Giele, Z. Kunzst and J. Winter. larXiv:0911. 19621 [hep-ph]. 

[38] R. K. Ellis, W. T. Giele, Z. Kunszt, K. Melnikov and G. Zanderighi, JHEP 0901, 012 (2009). 
[39] W.T. Giele and G. Zanderighi W. T. Giele and G. Zanderighi, JHEP 0806, 038 (2008). 
[40] A. Lazopoulos. larXiv:0911.5241 [hep-ph]. 
[41] A. Lazopoulos. larXiv:08123998 [hep-ph]. 

[42] P. P. Mastroha, G. Ossola, C. G. Papadopoulos and R. Pittau, JHEP 0806:030 (2008). 

[43] C.F. Berger, Z. Bern, L. Dixon, F. Febres-Cordero, D. Forde, H. Ita, D.A. Kosower, D. Maitre 
Phys. Rev. D 78, 036003 (2008). 

[44] Z. Bern and D. A. Kosower, Nucl. Phys. B 379 (1992) 451. 

[45] Z. Bern, A. De Freitas, L. J. Dixon and H. L. Wong, Phys. Rev. D 66 (2002) 085002. 

[46] Z. Kunszt, A. Signer and Z. Trocsanyi, Nucl. Phys. B 411, 397 (1994). 

[47] Z. Bern, L. J. Dixon and D. A. Kosower, Nucl. Phys. B 437, 259 (1995). 

Z. Bern, L. J. Dixon and D. A. Kosower, Nucl. Phys. B 513, 3 (1998). 

[49] Z. Bern, L. Dixon and D. Kosower, Ann. Rev. Nucl. Part. Sci. 46, 109 (1996). 

[50] V. Del Duca, L. J. Dixon and F. Maltoni, Nucl. Phys. B 571, 51 (2000). 

[51] S. Catani, S. Dittmaier and Z. Trocsanyi, Phys. Lett. B 500, 149 (2001). 

[52] F. A. Berends and W. Giele, Nucl. Phys. B 294 (1987) 700. 

[53] F. Maltoni, K. Paul, T. Stelzer and S. Willenbrock, Phys. Rev. D 67, 014026 (2003). 

[54] T. Stelzer and W. F. Long, Comput. Phys. Commun. 81, 357 (1994) [arXiv:hep-ph/9401258] . 

[55] S. Catani and M. H. Seymour, Nucl. Phys. B 485, 291 (1997) [Erratum-ibid. B 510, 503 
(1998)]. 

[56] S. Catani, S. Dittmaier, M. H. Seymour and Z. Trocsanyi, Nucl. Phys. B 627, 189 (2002). 
[57] Z. Nagy, Phys. Rev. D 68, 094002 (2003). 

28 



[58] J. M. Campbell and R.K. Elli s, Phys. Rev. D 62, 114012 (2000). The MCFM program is 
publicly available from http://mcfm.fnal.gov 

[59] J. M. Campbell, R. K. Ellis and F. Tramontano, Phys. Rev. D 70, 094012 (2004). 

[60] J. M. CampbeU and F. Tramontano, Nucl. Phys. B 726, 109 (2005). 

[61] S.D. Ellis and D.E. Soper, Phys. Rev. D 48, 3160 (1993). 

[62] M. Dasgupta, L. Magnea and G. P. Salam, JHEP 0802, 055 (2008) |arXiv:0712.30T4l [hep-ph]] . 

[63] T. Aaltonen et ai, The CDF collaboration public note 9850; See 
|http: //www-cdf ■fnal.gov/physics/new/top/2009/xsection/ttj_4 . linvf b/| 

[64] C. W. Bauer and B. O. Lange. larXrv:0905.4739l [hep-ph]. 

[65] C. F. Berger et ai, Phys. Rev. D 80, 074036 (2009). 

[66] K. Melnikov and G. Zanderighi, l arXiv:0910.367Tl 

[67] Contribution by S. Hoche, J. Huston, D. Maitre, J. Winter and G. Zanderighi to T. Binoth et 
at "The SM and NLO multileg working group: Summary report" , larXiv:1003.124l[ hep-ph]. 

[68] J. H. Kiihn and G. Rodrigo, Phys. Rev. D 59, 054017 (1999); Phys. Rev. Lett. 81, 49 (1998). 

[69] A. Alves, O. Eboh, T. Plehn and D. Rainwater. Phys. Rev. D 69, 075005 (2004). 



29 



